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It is well known that the ultra-relativistic Thomas-Fermi equation, amply adopted in the study 
of heavy nuclei, admits an exact solution for a constant proton distribution within a spherical core 
of radius Re- Here exact solutions of a generalized ultra-relativistic Thomas- Fermi equation are 
presented, assuming a Wood-Saxon-like proton distribution and its further generalizations. These 
solutions present an overcritical electric field close to their surface. The variation of the electric 
fields as a function of the generalized Wood-Saxon parameters are studied. 

PACS numbers: 97.60.Jd,04.70.-s,21.65.Mn,03.75.Ss,26.60.Dd,31.15.ht 

INTRODUCTION 

To study the electrodynamic properties of the bulk matter at nuclear densities a step proton distribution has been 
chosen Using the Migdal et. al. approximation, the ultra-relativistic Thomas-Fcrmi model, which governs 

this problem, reads 

^ = 0(.f-e(~a:), (1) 

where the proton density, Up, and the Coulomb potential at the center, V^(0), are given by 

x = k{r- R,), (2) 

k = 2V^(7r/6)i/6ny3, (3) 
6^(0) = (3^2„p)i/3. (4) 



The equation ([T]) admits the exact solution 



/ 1-3 [1 + 2-1/2 sinh(a-V3a;)] 

{ (x+b) ' ^ ^ " 

where integration constants a and b are: sinha = ll\/2, a — 3.439; b ~ (4/3)a/2-[3- 



(5) 



GENERALIZED ULTRA-RELATIVISTIC THOMAS-FERMI EQUATION 

In this section we want to look for exact solutions to a generalized ultra-relativistic Thomas-Fcrmi equation 



where 



[xY-fpOi-x), (6) 



fp{xb) — > 0, < < oo 

/p(-cx)) ^ 1, (7) 
f'p{x) < 0, for all X 

It is possible to write several distinct infinite 6-dependent sets of analytic solutions to the Thomas-Fermi Eq. 
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- Set 1 



{x-M = { V ^<^^' (8) 

^ ' ' ' for x>Xb, ' 



and leads to the following set of proton profiles 



fp{x; &) = I ^ ~ arctan(6a;))^ - ^(^g^.). , for x < xt, 



0, for X > Xb 



where a, (3 are real constants given by 



f3 = ^^^^±f^<l>{xb;b) - Xb, 



(10) 



because of the continuity of 4>{x; b), 4''{x; b) in Xb- 
The electric field for x < Xb,is given by 



2 2T^/^^2 1 ^ 



^(-;^)=(^eMo)^-TT(^- (11) 

The parameter 6 describes the width 25 (in cm) of the transition layer near the edge of the core 

Precisely 28 is the width of the transition layer of the core in which the electric field goes from its maximum to the 
half of its maximum. Now, let be be the value of b such that the electric field E{x = 0; be) is equal to the critical field 
Ec. Then 

bc;Sen ~ TT^TT-^, (13) 



and 

1 27 



229/6 5 



mc 



aony^icm). (14) 



where E„iax is the electric field at a; = to the step-proton distribution. We see that 5c can be of the order of the 
Bohr radius ao i.e. of order of 10^ electron Compton length. 



- Set 2 

^/^.M _ / i - |tanh(6a;), 

for X > Xb 



/3+x' 



and leads to the following set of proton profiles 

v3 



J- / j,\ _ / (^(1 — tanh(6a;)) — 6^ tanh(6a;)(l — (tanh(6a;))^), for x < Xb, 
jp[x;o)-<^^^ for x>Xb, ^ ' 



where a, (3 are real constants given by 

f _ J-^i tanh(bzj,) x 

J " — i,2\2 l+(tanh(6a:i,))2 J' 
anh(6xb) 
-(tanh(6x(,)) 



R — ( tanh(ba;fc) x _ 
P — U2-l-(tanh(5x!,))2 J ■^b, 



(17) 



because of the continuity of (t){x] b), <j)'{x; b) in Xb- 
The electric field for a; < a;6, is given by 



E{x; b) = 77r4T7^eV(0)4 . (18) 

^ ' (37r)i/2 ^ ' 2cosh^{bx) ^ ' 
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FIG. 1: The vaxiation of the potentiai, the proton density and the electric field strength near the edge of the core of the Set 1 
are plotted as functions of x, for a proton number density at the center rip = 0.57 • 10^^ cm~*. The dashed curve represents the 
function (j){x;bc); the dotted curve represents the distribution fp{x;bc); the solid curve represents the ratio E(x;bc)/Ec. 

The parameter b as above, describes the width 2S (in cm) of the transition layer near the edge of the core. Precisely 

0SET2 — OSETI ^ • (!») 

Now, let he be the value of h such that the electric field E{x = 0; be) is equal to the critical field E^. Then 

bc;Set2^^^, (20) 



and 



We note that 



0c;Set2 = 7^ Oc;Setl{cm). (21) 



^ - tanh(6a;) = ^775^ (22) 

2 2 1 + e(2''^) ^ ' 



which is well known in nuclear physics as Wood-Saxon profile. 
The Wood-Saxon profile can be generalized by 



(t>{x;a,b) = \ [i+ae(-^-)]/ (23) 

I for X > Xb, 
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with the following set of proton profiles 



fp{x;a,b)={ Y [l+aei-"-)]'^' j (l+ae-'.-)i/»(l+ae-'-)2' ^^^^ ^ ^ 

0, for X > Xb, 



(24) 



where a, /3 are real constants given by 



a 



(l+ae-'^^t)!/ 
- [(l+ae-^-.)^7°-l][l+ae-^-^ 



[(l+ae"''="'>)i/°-l][l+ae"''^i'] 
be-'"'b 



(25) 



because of the continuity of (t>{x] a, b), 4>' [x] a, b) in a;;,. 
We have 



hence the maximum of i?(a;; a, 6) is 



i?(x = 0; a, 5) = ^ , ^ . E,^a.. (27) 
(1 + a)^/'^(l + a) 



Set 3 

A(r..h\ - / 5- T2sinh"'(&a^), 
and leads to the following set of proton profiles 



0(:.;6) = ( V^^^"^"'^^")' (28) 



/(x;6) = < ('•2 72 '^inh (5a;)) 72(i+f,2'^2)3/2 , for x < Xb, ^2g) 
1 0, for X > Xb, 

where a, (3 are real constants. 

The Set 1, 5'ei 3 of analytic solutions to the Thomas- Fermi equation ^ belong to the more general following set 

, , / ci[$(x;a,6) +C2], for x < Xh, 
0(a:;a,6) = |_^^ for x > x„ ^^O) 



where $(x, a,6), Ci, C2 are given by 

(31) 



$(a;; a, 6) = ^i^i;2(l/2, a; 3/2, -&V), 

ci (lim.^^_fe7jc "J" (a;, a, fe) + lima.^oo(^' (a;, a, b))y^ , C2 = lim^^^oo ^ {x, a, 6) . 

and Fi;2 is the Gauss hyper-geometric function. For positive integer (> 2) or positive half-integer (> 3/2) values of 
Fi;2 can be written in terms of elementary functions (Table 
Also the Set 2 belongs to the more general set given by 

, , J ci[$(x;a,6) +C2], for x < Xb, 
0(a:;a,6)- (32) 

where 

1 $(x; a, 6) = j:,^^ A(i _ tanh(6y)2)"dy, ^ ^^^^ 

I ci = (limj;^_fcflc (2^, a, fc) + linij;^oo $ (a:^, a, 6)) \ C2 = lim^^^oo (a;, a, b) . 

For positive integer (> 1) or positive half-integer (> 1/2) values of a, can be written in terms of elementary 
fimctions (Table [Til. 

These results, obtained by explicit analytic formulae, complement the numerical results presented in 



* Electronic address: IrufEniQicra.itl 
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TABLE I: Generalized exact solutions to Set 1 and Set 3 



a 


$(x;a, 6) 


3/2 

2 
5/2 

3 
7/2 

4 
9/2 


— sinh~^(6a;) 
— ^ arctan(6a;) 

X 

(8(1+6^2^^)) 8f,^rctan(&a;) 

(15{l + 62i2)3/2) 

(24(1+6^2,^)^) (24(i+6!ia:2)) arctan(6x) 

(106(1+62 a!2)5/2) 



TABLE II: Gciicralizod exact sohitions to Set 2 



a 


$(a;;a, 6) 


1/2 


— sinh~"'^(tanh(&a::)) 


1 


-i tanli(6x-) - i ln(tanh(6a;) - 1) + 3 ln(tanh(6x) + 1) 


3/2 


— i tanli(6x)-^l — tanh(6a;)^ — | sinh~^(tanh(6a;)) 


2 


tanh(ba;)^ - \ tanli(6a;) - | ln(tanh(6a;) - 1) + | ln(tanh(6a::) + 1) 


5/2 


tanh(6a;)y^l - tanh(6a;)2 - ^ tanh(6a:)(l - tanh(6a:)^)^/^ - ^ sinli"^(tanh(6a:)) 
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